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Interaction metrics developed for the process control industry are used to evaluate 
decentralized control of sound radiation from bays on an aircraft fuselage. The metrics 
are applied to experimentally measured frequency response data from a model of an 
aircraft fuselage. The purpose is to understand how coupling between multiple bays of 
the fuselage can destabilize or limit the performance of a decentralized active noise control 
system. The metrics quantitatively verify observations from a previous experiment, in 
which decentralized controllers performed worse than centralized controllers. The metrics 
do not appear to be useful for explaining control spillover which was observed in a previous 
experiment. 


Introduction 

A CTIVE control of sound radiation due to tur- 
bulent boundary layer excitation of an aircraft 
fuselage has been the subject of much recent research. 
This problem is typically treated as one of minimizing 
sound radiation from individual bays of the aircraft 
fuselage, where a bay is a rectangular portion of the 
fuselage skin bounded by ring frames and stringers. 
The random, distributed nature of turbulent bound- 
ary layer (TBL) excitation means neighboring bays are 
driven by a spatially and temporally uncorrelated forc- 
ing function. 

The most common active control approaches assume 
feedback of a structural sensor such as an accelerome- 
ter or piezoceramic transducer . 1-4 Feedforward con- 
trol is less common in this application due to the 
absence of a reference signal which is correlated with 
the TBL excitation. Numerous analytical studies have 
demonstrated the potential for feedback control on this 
problem . 3,5-7 The feasibility of the approach has also 
been demonstrated in various experiments . 1,8-10 

Although active control is promising on individual 
fuselage bays, significant implementation issues must 
be addressed before a system can be deployed on an 
aircraft containing several hundred bays. For example, 
the control system components should be as simple 
and reliable as possible. Equally important is the 
question of a centralized versus a decentralized control 
scheme. For this application, the question is whether 
feedback control is implemented using only local infor- 
mation, or if a centralized controller coordinates mul- 
tiple feedback loops. A decentralized control system 
which relies only on local information is generally pre- 
ferred for reasons of simplicity and fault tolerance . 11 
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Much of the recent research in controlling sound ra- 
diation from panels has implicitly addressed this ques- 
tion by studying control systems which use matched, 
collocated actuator and sensor pairs . 2,3, 12 A matched, 
collocated actuator and sensor pair is physically collo- 
cated on the structure and has the same function form 
for the modal coupling coefficients . 13 This arrange- 
ment can be guaranteed to be dissipative, or passive, 
meaning it is inherently stable and tolerant to changes 
in the system being controlled. A matched, collocated 
controller is suitable for decentralized control because 
there is no need to account for destabilizing inputs 
from neighboring controllers. In spite of these ben- 
efits, this approach is categorized as a low- authority 
approach, since in practice passivity often severely lim- 
its the performance achievable through active control. 

Higher authority approaches based on linear 
quadratic optimal or robust control designs have also 
been studied for this problem . 1,8, 10 These approaches 
lack the robustness and stability guarantees of a 
matched, collocated controller, and are more compli- 
cated, but generally have better performance . 3 If de- 
centralized control is to be implemented with a higher 
authority approach, the design of the feedback loops 
may have to consider the destabilizing effect of inputs 
from neighboring controllers. For example, in a recent 
experiment comparing centralized and decentralized 
controllers on an aircraft-style sidewall, noise reduc- 
tion performance was limited by stability issues in the 
implementation of decentralized designs . 9 

The present paper considers the design and imple- 
mentation of decentralized controllers on an aircraft 
sidewall when the control system does not involve a 
matched actuator /sensor pair. Experimentally mea- 
sured frequency response functions are examined using 
interaction metrics which have been previously de- 
scribed in the process control literature . 11, 14 These 


l of 10 


American Institute of Aeronautics and Astronautics 



interaction metrics highlight frequencies where decen- 
tralized controllers could destabilize one another due 
to coupling in the system being controlled. The desta- 
bilizing effects could be mitigated either in the design 
of the individual decentralized control loops, or by 
implementing some degree of centralization, or coor- 
dination, among multiple feedback controllers. 

Because of the complex dynamics of a built-up struc- 
ture such as an aircraft fuselage, the analysis here uses 
data measured from a physical model of a fuselage side- 
wall. This model was constructed in order to study 
the control of sound radiation due to TBL excitation 
of an aircraft fuselage. This ensures the data has some 
of the complexities of an actual aircraft fuselage, such 
as varying boundary conditions due to manufactur- 
ing variabilities. But, it also means the analysis will 
include dynamics which may be particular to this spe- 
cific structure. Nonetheless, the analysis is useful for 
illustrating the interactions and control design consid- 
erations which may be found on a real fuselage. 

The paper begins with a description of a decentral- 
ized control system, and suggests two metrics for ana- 
lyzing interactions between feedback controllers. The 
experimental setup, which consists of a fuselage- style 
panel in a wind tunnel, is described next, including 
the data which were measured as part of the previous 
test. 9 The interaction metrics are then applied to this 
data and compared with the known panel dynamics 
in order to highlight areas where inter-bay coupling 
could present problems for a decentralized controller. 
The paper closes with a discussion of closed loop data 
for two control designs which had substantial spillover, 
possibly due to unmodeled interactions between neigh- 
boring feedback loops. 

Analysis 

This section introduces notation useful for analyzing 
the decentralized control problem and discusses two in- 
teraction metrics. Decentralized control is a common 
topic in controls textbooks, but interaction metrics de- 
scribed in these textbooks have not previously been 
applied to measured data from an aircraft-style struc- 
ture. This section summarizes the details needed to 
understand the metrics, and closely follows the excel- 
lent discussion of decentralized control given in Chap- 
ter 14 of Morari and Zafiriou. 14 

A block diagram of a decentralized control system 
is shown in Fig. 1. The G(z) block denotes a plant 
with m inputs and m outputs, where iq denotes the 
ith input and yi denotes the ith output. Feedback 
controllers are designated by the ki{z ) blocks, where 
the ith controller receives input from only the ith plant 
output and drives only the ith plant input, as shown 
in the figure. 



Fig. 1 Decentralized Control of plant G(z) 

Given this decentralized control topology, an im- 
portant question is the extent to which interactions 
between neighboring actuators and sensors will desta- 
bilize or reduce the performance of the control system. 
At the simplest level, this question concerns the mag- 
nitude of off-diagonal terms in the plant, G, relative to 
the diagonal terms. For example, if the transfer func- 
tion from the ith actuator to the ith sensor is denoted 
gu, then the off-diagonal term from the j th actuator to 
the ith sensor would be denoted gji . If the off-diagonal 
terms are large relative to the diagonal terms, it is 
easy to envision a scenario where decentralized con- 
trol loops could fight one another and eventually go 
unstable. 

A decentralized control scheme allows for lower or- 
der models, simpler control design techniques and has 
a significant advantage in implementation as there is 
no need to provide for signal paths between adjacent 
panel bays. Hence the goal is to design and imple- 
ment the individual feedback loops as if the plant G 
were diagonal. Interaction metrics can help under- 
stand how coupling, due to off-diagonal elements of 
G, constrains the design of the decentralized feed- 
back loops. 11,14 Significant coupling could constrain 
the loops to be very weak controllers in order to pre- 
serve stability of the fully coupled system. Although 
some insight into this question is available by examin- 
ing the cross terms of a frequency response function, 
this quickly becomes unwieldy for multiple feedback 
loops and does not yield direct design guidelines. 

Relative gain array metric 

This leads to the first interaction metric to be dis- 
cussed, which is the relative gain array (RGA). RGA 
can be used to quantify the diagonal dominance of a 
plant such as G in Fig. 1. The RGA is defined as 14 

A(G(w)) = G(w) x (G(u)~ 1 ) t (1) 

where x denotes the element by element product. 
Thus the RGA of a frequency dependent square matrix 
G is the same size as G and is also frequency depen- 
dent. A few properties of the RGA which are relevant 
for analyzing a decentralized controller are: 

• A(G) is independent of scaling on the inputs and 
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outputs of G 


G-G 


• A (G) is close to the identity matrix for plants 
which are diagonally dominant. The RGA num- 
ber, or sum norm of the RGA matrix, quantifies 
the closeness of the RGA to the identity matrix : 11 

RGA number = ||A - I\\ 8um ( 2 ) 

• if G is stable and each decentralized controller is 
stable, then A (G) = / implies stability of the en- 
tire system . 11 



When the plant is diagonally dominant, it is unlikely 
that the operation of one decentralized controller will 
have any impact on the operation of any other con- 
troller, which is obviously desirable. The RGA has a 
simple form for a (2 x 2) system: 


A = 


An 

1 — An 


1 — An 
An 


( 3 ) 


where A n =^11^22/(^11^22 -^12^21). 

Note that the RGA can be computed from a fre- 
quency response function description of a plant. This 
means some insight into the diagonal dominance of a 
plant can be obtained simply by examining measured 
input/output data. There is no need to compute a re- 
alization (e.g. a state-space or ARMA model) of the 
plant under consideration. Because the RGA is scaling 
independent, a simple test for diagonal dominance is 
to look for elements of A (G) which are much greater 
than one, or equivalently where A (G) is significantly 
different from the identity matrix. 

The RGA is most commonly discussed in the context 
of the pairing problem, where one must decide which 
inputs to pair with which outputs . 11, 14 It is useful 
in that context because it only needs to be computed 
once for the candidate input /output set, and different 
actuator /sensor pairings can be studied without hav- 
ing to recompute the RGA . 11 Here, however, the desire 
to have independent control loops on each fuselage bay 
uniquely determines the actuator /sensor pairing. 

Although the RGA is useful for analyzing the diag- 
onal dominance of a plant, it does not provide bounds 
on the allowable control gains. 

fi metric 

A more detailed interaction metric which does pro- 
vide gain bounds can be derived using concepts from 
robust control analysis. Following the approach de- 
scribed in Morari and Zafiriou , 14 the block diagram 
of the decentralized control system can be redrawn as 
shown in Fig. 2 . This diagram is equivalent to the 
first diagram, but it illustrates how the off-diagonal 
elements of G can be viewed as perturbations on the 


Fig. 2 Decentralized control schematic for inter- 
action analysis 

diagonal elements. The control loops, ki , are designed 
for the diagonal plant, G. The stability of the resulting 
decentralized controller, IT, on the full system, G, can 
then be quantified by the robustness of the controller 
with respect to the additive uncertainty, (G — G). The 
additional signal, d, in the figure corresponds to a dis- 
turbance input measured by the sensors, and n denotes 
random sensor noise. 

At this point it is helpful to introduce a few stan- 
dard quantities for analyzing the control system. The 
sensitivity of the system in the figure is defined as a 
transfer function from d to 7 /, 15 and is written 

E= T^gk 


Similarly, the complementary sensitivity, which is the 
transfer function from n to ?/, is written 


H = 


GK 

I -GK 


where the relationship 


( 5 ) 


E + H = / 


(6) 


holds . 15 These quantities have been defined assuming 
positive feedback in order to be consistent with the 
experimental results discussed later in the paper. 

Some general, well-known 15 properties of these 
transfer functions are worth mentioning. For good 
suppression of the disturbance, d, the sensitivity E, 
should be small where d is significant. For a MIMO 
control system, this requirement can be written in 
terms of the maximum singular value of E , denoted 
ct(E), as 

*(£)«! (?) 

Similarly, to minimize the effect of the noise n on the 
response y , the complementary sensitivity H should 
be small: d(H) <C 1 . 
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The return difference for the system in Fig. 2 can be 
factored as 11,14 


(/ - GK) = (I- L h H ) (/ - GK) (8) 

where 

L h = (G — G)G~ 1 (9) 

quantifies the relative error in approximating G by the 
diagonal G. The second term on the right hand side 
of Eq. 8 is the return difference if the plant is assumed 
to be diagonal (i.e. G = G). The complementary 
sensitivity for the diagonal plant is written 


I -GK 


(10) 


which is itself diagonal. The ith diagonal of H can be 
written 


Qii^i 

1 


(ii) 


The full H is stable if each loop, hi , is stable. If H is 
stable, then the full control system is stable if and only 
if the Nyquist criterion is satisfied for the first term on 
the right hand side of Eq. 8: (/ — LhH). 14 

From the small gain theorem, a sufficient condition 
for stability is to limit the spectral radius of LhH , 
denoted p(LnH ), such that 


p(L h H) < 1 (12) 


This is satisfied if the product of induced norms is less 
than unity, 

\\Lh\\\\H\\ <1 ( 13 ) 

This equation provides a simple measure for design 
purposes. Lh is strictly a function of the open loop 
plant so it is independent of controller design. Fur- 
thermore Lh can be computed from experimental 
frequency response data without requiring a system 
realization. If H is diagonal, \\H\\ = max|/^|, pro- 
vides a gain bound for the controller that need only 
consider the most aggressive loop. 

However, the simplicity of the metric in Eq. 13 
comes at a price of conservativeness. The norm con- 
dition in Eq. 13 is only sufficient for the Nyquist cri- 
terion, and therefore can be quite conservative. One 
way to reduce this conservatism is by applying a fre- 
quency dependent scaling to the norm calculation on 
the normalized plant error, Lh • This yields a struc- 
tured singular value, 16,17 or //, condition for stability 


\\H(lu)\\ <im-\L h (lj)), 0<uj<n/T (14) 

where //(•) yields a real, positive scalar, p is computed 
for a particular uncertainty structure, which for decen- 
tralized control corresponds to the diagonal structure 
of H. The use of Eq. 14 in place of Eq. 13 can reduce 
conservatism considerably. 
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Fig. 3 Test panel schematic 


Experimental Setup 

The interaction metrics were applied to data col- 
lected in a previous controls experiment. 9 This inter- 
action analysis was done by re-examining frequency 
response and closed loop data from the earlier experi- 
ment. 


Test Panel 

The test configuration has been described in detail 
elsewhere 1,9 and is summarized here. The test struc- 
ture consists of an aluminum panel which models the 
sidewall of an aircraft fuselage, absent the fuselage cur- 
vature. A schematic of the panel is shown in Fig. 3. 
The structure consists of six bays separated by vertical 
aluminum frames and horizontal stringers. The dimen- 
sions of each bay are 50.8 cmx 25.4 cmx 0.16 cm. The 
bays are referred to in later sections of the paper by 
the numbers indicated on the schematic. In-plane ten- 
sion was applied to the panel in the vertical direction 
to simulate the hoop stress in a pressurized aircraft 
flying at 12,000 m. 

The panel was mounted in the sidewall of the Struc- 
tural Acoustic Flow Apparatus (SAFA) at NASA Lan- 
gley Research Center. SAFA is a low speed wind 
tunnel which is used to study turbulent boundary layer 
excitation of panels. A photo of the panel installed in 
the tunnel sidewall is shown in Fig. 4. 

Transducers 

The control system on each bay consisted of a sin- 
gle actuator and four accelerometers whose outputs 
were summed and treated as a single sensor. The ac- 
tuator was a center- mounted piezoceramic (PZT-5A, 
15.2 cmx 10.2 cmx 0.038 cm). The sensor response 
on each bay was created by summing the response 
of four accelerometers which were arranged in a di- 
amond pattern centered on the PZT. The four ac- 
celerometers in the diamond array on each bay are 
indicated in the schematic by blackened circles. The 
output of this sensing configuration is an approximate 
measure of the response of the first radiation mode 
of the panel. 1 Eleven additional accelerometers were 
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Fig. 4 Test panel mounted in wind tunnel 


mounted on each bay and used in a separate data ac- 
quisition system to estimate the radiated sound power 
from each bay. 

The decentralized control topology considered here 
consisted of independent single input /single output 
(SISO) control loops running on each bay. The cen- 
tralized control topology consisted of a fully cou- 
pled, multiple-input/multiple-output (MIMO) topol- 
ogy, where the number of inputs and outputs equaled 
the number of bays being controlled. 

Frequency response data between actuators and sen- 
sors were collected using the digital control system. 
This ensured an accurate measurement of all gains and 
delays in the control path. The sample rate for all data 
was 3000 Hz. 


Results and Discussion 

The interaction metrics were used to analyze decen- 
tralized control of two pairs of bays on the test panel. 
The first pair consisted of bays one and two, and the 
second pair consisted of bays one and four. A more 
realistic interaction analysis will have to include more 
than just two bays at a time, especially if decentralized 
control is to be implemented on an aircraft fuselage. 
However, only two bays at a time are considered here 
to simplify interpretation of the results. 

The full plant matrix, G, for each two-bay configu- 
ration was either 


9 li 9 12 
921 922 


or 


9 li 9 14 

94,1 94 , 4 , 


(15) 


depending on which pair of bays was being studied. 
One obvious centralized controller for this system is a 
two input/two output (or MIMO) configuration which 
considers off-diagonal terms of G in the control design. 
A decentralized configuration consists of two single 
input /single output (SISO) controllers which neglect 


off-diagonal terms in the design of the individual SISO 
control loops. 

Experimentally measured values of the matrices in 
Eq. 15 as a function of frequency are shown in Fig. 5. 
The figures show measured frequency response func- 
tions (FRFs) in terms of average acceleration, across 
the diamond array, per unit voltage input to the PZT. 
The data on bays 1 and 2, in Figs. 5(a) and 5(b), were 
collected by simultaneously driving the PZTs on those 
bays with independent random noise sequences. Simi- 
larly, the data on bays 1 and 4, in Figs. 5(c) and 5(d), 
were collected by simultaneously driving the PZTs on 
bays 1 and 4 only. In each plot, the solid line shows 
the FRF from the PZT to the sensor on the same bay, 
while the dashed line shows the sensor response on the 
neighboring bay. The frequency range in the plots is 
restricted to 150 Hz to 700 Hz for simplicity. Coupling 
occurs outside this frequency range, but the selected 
range is of particular interest for active control. 

A simple check for diagonal dominance of each (2x2) 
plant is to examine the separation between the solid 
and dashed lines in Fig. 5. Frequencies where the two 
lines coincide, or where the dashed line lies above the 
solid line, indicate regions where the plant is not di- 
agonally dominant. Several such areas are visible in 
the plots. For example, the solid and dashed lines co- 
incide near 180 Hz in all four plots. For bays 1 and 2, 
Figs. 5(a) and 5(b) show coupling near 325 Hz , 425 Hz, 
and 625 Hz. On bays 1 and 4, coupling is evident near 
260 Hz, 325 Hz, 375 Hz, and 425 Hz. 

The operating deflection shapes offer some insight 
into the panel dynamics at these frequencies. The de- 
flection shapes at four frequencies of interest, due to 
a random input to the PZT on bay one, are shown 
in Fig. 6. These were obtained by plotting the real 
part of the phased responses of the fifteen accelerome- 
ters on the panel. They reveal that many of the FRF 
peaks correspond to resonances of multiple, coupled 
bays. For example, 180 Hz, in Fig. 6(a), corresponds 
to be a (3,1) entire panel mode, with a nodal line on 
each of the vertical panel frames. Because it is a panel 
mode, there is significant coupling between all bays 
on the panel. At 260 Hz, the driven bay and its verti- 
cal neighbor each vibrate in their (1,1) bay modes, but 
180° out of phase with one another. There is consider- 
ably less coupling between the side by side bays at this 
frequency. For example, the response of bay 2 when 
bay 1 is driven (Fig. 5(a)) is considerably lower than 
the response of bay 4 when bay 1 is driven (Fig. 5(c)). 

Coupling between side by side bays is present above 
600 Hz. The peak in Figs. 5(a) and 5(b) at 625 Hz cor- 
responds to the resonance of the (5,1) bay mode of bay 
one. The peak at 641 Hz in those figures is the (5,1) 
bay mode of bay two. The FRF data shows significant 
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Fig. 5 Frequency responses from measured data: (a), (b) bays 1 and 2; (c), (d) bays 1 and 4 


coupling between these bays in this frequency range. 

Although the FRFs are informative, the interaction 
metrics offer more concise information on bay inter- 
actions. The RGA number from Eq. 2 is shown in 
Fig. 7 for bays one and two. Wherever the RGA 
number is non-zero, the RGA matrix differs from the 
ideal identity matrix and the plant G is not diagonally 
dominant. The plot shows that for much of the fre- 
quency range, the diagonal terms of G are dominant. 
The only areas where the number deviates significantly 
from zero are near the panel mode at 180 Hz and 
above 600 Hz. The coupling near 180 Hz could be 
problematic for active control of sound radiation. Low 
order panel modes, such as at 180 Hz , are more effi- 
cient sound radiators than higher order panel modes. 
If the authority of the decentralized controllers have 
to be reduced here to maintain stability, noise reduc- 


tion performance could be impacted. The coupling 
near 600 Hz may be less of a problem because it is 
associated with a higher order mode. The gain of the 
controllers could be rolled of above 600 Hz to avoid 
the possible destabilizing affects of this coupling. 

The RGA number for bays one and four is shown 
in Fig. 8. The data indicates the plant matrix is not 
diagonally dominant below around 450 Hz. In addi- 
tion to the peak near 180 Hz , there are several large 
peaks between 240 and 275 Hz. These peaks generally 
agree with the FRF data, which suggested significant 
coupling between these bays in this frequency range. 
This coupling between the (1,1) bay modes at 260 Hz 
can also present a problem for decentralized control 
of sound radiation, as with the 180 Hz mode. The 
RGA number is close to zero above 600 Hz , which is 
expected since the FRFs indicated very little two way 
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Fig. 6 Operating deflection shapes due to input on bay one. 


coupling between these bays in this frequency range. 

To summarize the RGA data for these bay pairs: the 
plant (G) for bays one and two is diagonally dominant 
across a large frequency range except near the panel 
mode at 180 Hz and above 600 Hz. In contrast, the 
plant for bays one and four is diagonally dominant only 
for frequencies above about 450 Hz. If decentralized 
control is to be used on either bay pair, the possible 
destabilizing effects of these areas where the plants are 
not diagonally dominant should be considered. This 
appears to be most important for bays one and four. 

The interaction metric fi^Ln) is plotted for both 
bay pairs in Fig. 9. The solid line corresponds to bays 
1 and 2, while the dashed line corresponds to bays 1 
and 4. As shown in Eq. 14, h(Lh) provides a numeric 
constraint on the diagonal complementary sensitivity, 
H. 

For bays one and two the /i data generally agrees 
with the RGA data. /r(L#) for these bays is greater 
than unity near 180 Hz and above 620 Hz. For bays one 
and four, h(Lh) indicates some additional frequency 
regions where the plant may not be diagonally domi- 
nant. For example, /r(L#) peaks near 225 Hz but only 
slightly exceeds unity at 260 Hz. In contrast, the RGA 
peaked at 260 Hz , and was unity at 225 Hz. fi( L h ) also 
exceeds unity at 320 Hz and 380 Hz , and then peaks 
at 420 Hz. 

To summarize, the /r(L#) metric confirms that sta- 
ble decentralized control of bays one and two can be 
achieved with relatively little constraint on H. For 
bays one and four, however, fi(Ln) indicates several 
frequency ranges where H is constrained, if stability 
of the decentralized controller is to be guaranteed. 



Fig. 7 RGA number for bays 1 and 2 


Analysis of Closed Loop Data 

Two decentralized control designs from an earlier 
experiment are reexamined here. The first controller 
was designed for bays one and two on the panel, and 
the second controller was designed for bays one and 
four. In both control designs the off-diagonal terms 
of the plant were ignored. These particular controllers 
were chosen because they were stable but had substan- 
tial spillover of control energy. The closed loop results 
and the metric in Eq. 14 are compared in order to 
evaluate the metrics for this application. 

Before discussing the controllers in detail, typical re- 
sults for centralized and decentralized control of these 
bay pairs are summarized here. The results are listed 
in Table 1 in terms of the change in radiated sound 
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Fig. 8 RGA number for bays 1 and 4 



Fig. 9 h(Lh) for bays 1 and 2 (solid line); bays 1 
and 4 (dashed line) 

power with control relative to the no control case, on 
each bay of interest. The reductions were computed 
from the fifteen accelerometer responses, and were 
integrated from 150 to 800 Hz. The corresponding 
tunnel flow speed was Mach 0.125. The first column 
specifies the control type, either centralized (MIMO) 
or decentralized (SISO). The next two columns show 
sound power reductions measured from bays 1 and 2 
when those two bays were controlled. The last two 
columns show sound power reductions when only bays 
1 and 4 were controlled. The results demonstrate that 
for both bay pairs, the modest reductions obtained 
with a centralized design could not be obtained with 
the decentralized designs. As discussed by Gibbs, 8 9 
the decentralized designs would go unstable before 
achieving the reduction performance of the centralized 
designs. 

Both centralized and decentralized control designs 



Bay Number 

Control 

1 

2 

1 

4 

MIMO 

-3.0 dB 

-5.4 dB 

-3.0 dB 

-2.7 dB 

SISO 

-2.2 dB 

-3.1 dB 

-1.4 dB 

-1.4 dB 


Table 1 MIMO vs. SISO control results 9 


were computed using a predictive control formulation, 
as described previously. 9 The controllers were com- 
puted to minimize a cost function of the form 

J = EII^)II 2 + a|Ka : )|| 2 (16) 

k = 0 

The parameter p is the cost horizon over which the 
control solution is optimal. 18, 19 The solution approx- 
imates a linear quadratic regulator for large values of 
p. A state estimator is required when implementing 
the controller, so the resulting control implementation 
is very similar to an LQG control solution. 15 No fre- 
quency shaping of either the response or control effort 
terms in Eq. 16 was used. For the SISO (decentralized) 
control designs, this means that possible destabilizing 
effects of coupling were neglected. 

A high order state space model of the plant was 
computed using an observer Kalman filter realization 
algorithm. 20 For both MIMO and SISO control de- 
signs, the control effort parameter, A, in Eq. 16 was 
reduced to the smallest value that still produced a sta- 
ble controller when implemented. This approach may 
not produce controllers which are robust to plant vari- 
ations, but it provided a consistent design goal for the 
different control configurations. 

Auto-spectra of the sensor response on bay one with 
and without control are shown in Fig. 10 for decen- 
tralized control of bays one and two. The sensor is 
the summed response of the diamond accelerometer 
array. The solid line shows the auto-spectrum without 
control, and the dashed line shows the response with 
control. The open loop response is dominated by the 
panel mode at 180 Hz , but this was reduced by 10 dB 
by the control system. Spillover is evident at 300 Hz, 
350 to 380 Hz, and between 400 and 500 Hz. 

Although the RGA and the /jl(Lh) metric indi- 
cated potential coupling problems at 180 Hz and above 
600 Hz, the control system performed fine at these fre- 
quencies. To understand this, the more detailed metric 
in Eq. 14 was computed and is shown in Fig. 11 for this 
decentralized control configuration. The product of 
&(H) and /jl(Lh) should be less than unity to guaran- 
tee stability of the decentralized controller. The data 
in the figure show that this metric was satisfied except 
near 450 Hz and near 640 Hz. Thus even though the 
control design did not consider the constraint imposed 
by p(Lh), the resulting design happened to satisfy the 
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Fig. 10 Auto-spectrum of diamond accelerometer 
array on bay 1: Open loop (solid line); decentral- 
ized control of bays 1 and 2 (dashed line). 



Fig. 11 (j(H)/jl(Lh) for decentralized control of bays 
1 and 2 

constraint in Eq. 14 at most frequencies. Although the 
constraint violation at 450 Hz may explain the spillover 
seen at that frequency, no such spillover occurred at 
620 Hz , where the constraint was also violated. Sim- 
ilarly, the constraint in Fig. 11 does not offer insight 
into the spillover at 300 Hz , since it does not indicate 
a problem in that area. 

Results for the second control configuration are 
shown in Fig. 12. The figure shows the auto-spectra 
of the sensor response on bay one for decentralized 
control of bays 1 and 4. As with bays one and two, 
the controller produced good attenuation near 180 Hz 
but spilled over near 300 Hz. g(H)h(Lh) for this case 
is shown in Fig. 13. This constraint exceeds unity at 
more frequencies than for bays one and two. There is 
a peak in Fig. 13 at 300 Hz which could explain the 
spillover seen there. However, spillover is not seen in 


Fig. 12 Auto-spectrum of diamond accelerometer 
array on bay 1: Open loop (solid line); decentral- 
ized control of bays 1 and 4 (dashed line). 



Fig. 13 (j(H)h(Lh ) for decentralized control of bays 
1 and 4 

other frequency ranges where the constraint was vio- 
lated, such as at 225 Hz or 420 Hz. 

Conclusions 

Interaction metrics which were previously described 
for analyzing decentralized controllers 14 were applied 
to experimentally measured data from a model of an 
aircraft sidewall. Frequency response data from pairs 
of bays on a fuselage panel were examined. The Rel- 
ative Gain Array, or RGA, was studied as a means 
to quantify the diagonal dominance of a plant. The 
second interaction metric studied was based on a fi- 
analysis of the off diagonal terms of the plant relative 
to the diagonal terms. The /i metric provided a quan- 
titative constraint on the maximum singular value of 
the complementary sensitivity of a decentralized con- 
troller in order to guarantee stability. 
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On the panel studied here, the RGA and the /i met- 
ric showed that bays sharing a short vertical boundary 
were less coupled than bays sharing a long horizontal 
boundary. Both metrics also indicated that most cou- 
pling occurred when the resonance of one bay coupled 
into the resonance of a neighboring bay, as would be 
expected. For bays sharing a long horizontal bound- 
ary, satisfying the ft constraint in order to guarantee 
stability of the decentralized controller could greatly 
complicate the control design procedure. Further 
analysis is needed to accurately correlate inter-bay 
coupling to structural characteristics such as shared 
boundaries. 

Control designs from a previous experiment were 
found to violate the ft interaction constraint at sev- 
eral frequency points, although the controllers were 
stable when implemented. This verifies the conserva- 
tiveness of this constraint. Spillover produced by the 
controllers was not easily correlated with violations of 
the ft constraint. 

More work needs to be done on applying these met- 
rics to this problem. The next step is to use the ft 
constraint in the design of decentralized controllers. 
It will be important to determine if the resulting con- 
trol designs are too conservative to produce acceptable 
noise reduction performance. 
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